Variational solution of the T-matrix integral equation 
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We present a variational solution of the T-matrix integral equation within a local approximation. 
This solution provides a simple form for the T matrix similar to Hubbard models but with the local 
interaction depending on momentum and frequency. By examining the ladder diagrams for irre- 
ducible polarizability, a connection between this interaction and the local-field factor is established. 
Based on the obtained solution, a form for the T-matrix contribution to the electron self-energy in 
addition to the GW term is proposed. In the case of the electron-hole multiple scattering, this form 
allows one to avoid double counting. 

PACS numbers: 71.10.-w 



I. INTRODUCTION 

As a result of the first cycle of an iterative solution of 
the Hedin equations?i the commonly used GW approx- 
imation (GWA) models the electron self-energy as the 
product S = iGqWq of a noninteracting Green function 
Go and a dynamically screened Coulomb interaction Wq 
obtained within the random phase approximation (RPA). 
The GWA that describes the long-range screening well 
has been successfully applied to a broad spectrum of ma- 
terials where the interaction is not too strong and screen- 
ing effects dominate. However, the GWA encounters dif- 
ficulties (first of all in its description of the satellite struc- 
ture) in the case of systems with localized states where 
short-range interaction prevails For such systems, one 
has to use a theory beyond the GWA. This theory can be 
based on both an improvement of the RPA to get a more 
realistic screening picture and an inclusion into calcula- 
tions of the electron self-energy of the higher-order terms 
in the screened interaction. 

The first attempt to improve the RPA by including 
the effects of the exchange-correlation (XC) hole is well 
known to have been undertaken by Hubbard,'* who intro- 
duced the so-called local-field factor. The concept of the 
latter is that all corrections to the RPA can be formally 
reduced to it. However, the Hubbard local-field factor 
Q{c\} includes the frequency-independent exchange hole 
correction only. Diagrammatically such t/(q) can be ex- 
actly derived by summing the ladder diagrams for irre- 
ducible polarizability with a contact interaction and non- 
interacting Green functions (see, e.g., Ref.i). In order 
to explicitly include into consideration the full static XC 
hole around the screening electron, Singwi et al.^ have 
obtained more sophisticated expression for ^(q) which 
contains the equilibrium static pair-correlation function.^ 
Further essential improvements in the derivation of the 
local-field factor have recently been done by different au- 
thors (see, e.g., Refs. Iala il0|) who have studied the fre- 
quency dependence of the XC hole. 



The concept of the local-field factor has taken on a new 
physical meaning in time-dependent density-functional 
theory (TDDFT).*! In the TDDFT within hnear re- 
sponse theory, the dynamical factor C/(q, cj) is linked to 
the XC kernel fxc{f\,^)- The latter plays the role of the 
time-dependent (TD) XC interaction in addition to the 
Coulomb repulsion Vc- As a result, the response function 
R can be written as^^ 

R{q)^P{q)+P{q)v,{\q\)R{q), (1) 

where the irreducible polarizability P is defined by the 
equation 

P{q)^P%q)+P''{q)fxc{q)P{q). (2) 

Here and in the following we use the four-momentum 
variable g as a shorthand for (q, w). In Eq. Q P^ 
is the RPA irreducible polarizability and fxcil) = 
-'Vc{H)G{q). 

In order to derive Eq. jSJl from the Hedin equation for 
the irreducible polarizability^ 

^^'^^ = " (1^ / ^^^(^)G(^ - 9)' (3) 

where G{k) is the Green function and A(fc, q) is the 
vertex functionjii the latter must depend on one four- 
momentum q only (see, e.g., Refs. [T3 and ITsl) . i.e., 

A(fc, q) = -. . ! Noo^ ^ ' (4) 

1 - fxc{q)P^{q) 

which finally leads to P{q) = P"(g)A(q). 

Diagrammatically such a form for the vertex function 
has been obtained by Richardson and Ashcroft in Ref. ^ 
using a local approximation-'^^ within a variational ap- 
proach. They have summed an infinite number of self- 
energy, exchange, and fluctuation terms in the diagram- 
matic expansion of A. In contrast to the Hubbard t/(q), 
the local-field factor derived by this summation is a dy- 
namical one. 
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The representation of the vertex function allows 
one to include vertex corrections into t he calcu lation of 
the electron self-energy (see, e.g., Refs.|llli3). Thus, 
the concept of the local-field factor suggested by Hubbard 
considerably simplifies a problem of vertex corrections 
calculations in numerical applications and transfers all 
weight of the problem to calculations of the local-field 
factor (or XC kernel) for real systems. 

Fundamentally distinct way to go beyond the GWA is 
based on the use of the T matrixiiSiiS The T-matrix ap- 
proximation (TMA) originally was established to study 
strongly correlated fermion systems with short-range in- 
teraction and is strictly valid in the limit of an al- 
most filled or, because of particle-hole symmetry, an al- 
most empty band.-^ -^ This approximation allows one 
to include processes involving multiple scattering be- 
tween two electrons or two holes. This fact makes 
the TMA capable of describing a satellite structure, 
for example, in Nii^2i22iSi24 However, these calculations 
were performed using either a statically screened model 
interaction^ or the Hubbard parameter U within Hub- 
bard models iSS42iiS242^ In the latter, the T matrix in mo- 
mentum space depends only on one four-momentum [as 
well as the vertex function Q expressed in terms of the 
local-field factor] and schematically can be represented 



as 



T{q) 



U 



1 - UK{q) 



(5) 



where K{q) is the Fourier transform of the product of 
two Green functions. In contrast to Eq. Q), an object of 
principal concern here is the local interaction U. 

Heuristically combining the simplification of Hubbard 
models, the T-matrix formalism of Ref . and a contact 
interaction W = >V(r,r';w = 0)i5(r - r') as in Ref.lH a 
GW + T matrix approach has recently been developed in 
Ref. This approach has effectively been applied to an 
excited electron lifetime in ferromagnetic Fe and Ni. In 
fact, comparing with the Hubbard models, one can find 
that the model short-range interaction U in the method 
of Ref. is replaced by the statically screened Coulomb 
interaction Wo(q, a; = 0). The possibility of such replace- 
ment was recently suggested by several authors 
Additionally, the importance of frequency dependence of 
the Hubbard U has been demonstrated in Ref. l29l 

The motivation of this work is to find a way that allows 
us to get the same result as the Hubbard model simpli- 
fication for the T matrix which is free of model parame- 
ters and with the momentum- and frequency-dependent 
local screened interaction. In order to accomplish this, 
we employ a variational method*'^° to solve the Bethe- 
Salpeter equation for the T matrix within a local approx- 
imation. As a result, the T matrix depends only on one 
four-dimensional wave vector, such as the vertex function 
expressed in terms of the local-field factor. 

The paper is organized as follows. In Sec^I we con- 
struct variational functionals and obtain from the van- 
ishing of their variational derivative a solution of the T- 
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FIG. 1: Feynman diagrams for T^^; (a) and T^^f" (b) in co- 
ordinate space. The T matrix is shown by the shaded square. 
The wiggly lines signify the dynamically screened Coulomb in- 
teraction W. The solid lines with arrows represent the Green 
function G. 



matrix integral equation. In order to connect this solu- 
tion with the results known from the literature, in Sec. lIIII 
we sum the exchange terms in the diagrammatic expan- 
sion of the irreducible polarizability by using the T ma- 
trix obtained. In Sec. IIVI we derive basic formulas for 
the electron self-energy beyond the GWA. Finally, the 
conclusions are given in Sec.lVl 



II. T MATRIX 

In this section we present mathematical expressions 
which lead to a simple form for the T matrix depending 
on a four-momentum only. We start from the T matrix as 
an object which will help us in our treatment of the ladder 
diagrams both for the irreducible polarizability P and 
for the electron self-energy S. The matrix is defined by 
the following Bethe-Salpeter equatiorii2iSii2& (Feynman 
diagrams are shown in Fig. 

T^,, (1, 2|3, 4) = W{1, 2)S{1 - 3)S{2 - 4) 

+ W{1,2) J dl'd2'K:,,{l,2\l',2') 

xr,",,(l',2'|3,4), (6) 

where W is the dynamically screened Coulomb interac- 
tion and a labels the spin, a can be specified as e — e 
in the case of multiple scattering between two electrons 
or holes and as e — /i in the case of multiple scattering 
between an electron and a hole. The kernel K^^, is the 
product of the Green functions 6*0^(1, 2): 

x:7(l,2|l',2')=iG.(l,l')G.K2,2'), 
X^-''(l,2|l',2') = zG.(l,l')G..(2',2). 

We have used the shorthand notation 1 = (ri,ti). As 
in the majority of practical schemes (including the com- 
monly used local-density approximation schemes), we 



3 




the right-hand part of the bubble K"^,) and e — h (down- 
FIG. 2: Feynman diagrams for T^^'f (a) and F^",*^ (b) in mo- directed arrow) cases, 
mentum space. 



suggest for simplicity that the system considered has 
properties of a homogeneous system. As a result, the 
T matrix © in momentum space has the formi^ 

X S[pi±p2-{P3±P4)]- (7) 

In the notations, we use the upper sign for the e — e 
and the lower sign for the e — h case. The (5-function in 
Eq. {Tj) reflects the conservation of total four-momentum 
in a homogeneous system and 

+ J dkWik)G^{piTk)G,'{P2 + k) 

^Ka'{PlTk,P2 + k\p3,P4)- (8) 

Feynman diagrams for F^^, are shown in Fig. |21 It is con- 
venient to introduce the total center-of-mass wave vector 
and the relative wave vectorsi^ 

Q =Pi±P2=P3±P4, ^{pi TP2), q' = ^(P3 TP4)- 

In terms of these new variables the function F"„, from 
Eq. ijHJl can be cast into the form given by 

KAq.'iM) 

= O (^Q + 9, T + 9', ±\q t q)- (9) 



Defining 



-G„{QTk)G,,{k) (10) 



and 



$^,,(q,fc,g) = 5{q-k) 



W{±qTk)KZ„,^Q{±-QTk), (11) 



one derives from the starting equation ||SJ) the relation 
j dk<^Za- {q, k, Q)K,, (fc, q',Q) = Wi±q T ?')• (12) 

The integral Eq. (|12|l can also be obtained from the 
vanishing of a functional derivative 



SF°'[G,W,r] 
SKAl^q\Q) 



0, 



(13) 



where F, a functional of three independent variables 
G, W, and F, is given by 

F"[G, VF,F] 

= dkdq'dQT"^^,{k,q',Q)K':^,^Q{±^QTk) 

X { J dp$^,,(fc,p,Q)F^^,(p,g',Q) 

-2Wi±kTq')}<.,,Qi±lQTq')- (W) 

Taking a trial solution in the spirit of the local approxi- 
mation of Ref. Is 



F^,,(g,q',Q)=0(Q), 



we find that 



(15) 



(16) 



where 



w:aq) - [K,,{Q)r'M:AQ)[K.'{Q)r\ 

M^,,{Q)^ I dqdpnZ,,Q{q)W{q~p)nZ,,n{p). 



Thus, we have obtained the T matrix as a function of 
the total center-of-mass wave vector Q only. Comparing 
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Eq. (|16|l with Eq. ((SJ , one can see that instead of the Hub- 
bard parameter U we have a momentum- and frequency- 
dependent local interaction W^^,{Q). The structure of 
F"^, in terms of this local interaction is schematically 
illustrated in Fig. 13 

III. IRREDUCIBLE POLARIZABILITY 

We will show here that the T matrix produces 
the irreducible polarizability in the form of Eq. with 
the local-field factor existing in the literature. Actually, 
the T matrix allows one to sum the all-order exchange 
diagrams in the irreducible polarizability diagrammatic 
expansion (corresponding Feynman diagrams are shown 
in Fig. 0: 

P(l,2)= P°(l,2) 

+ 51 / f^3d4d5d6G^(l,3)Ga(4, 1) 

X T,V"(3,4|5,6)G,(2,6)G,(5,2). (17) 
In momentum space, we have 

Pip) = - E / dkdq^l-^ik){Sik - q) 

+ K-'\k + ^p,q+^p,p)Kl-'^^{q)}. (18) 

By substituting the T matrix Hlt)|) into Eq. H18|l one 
obtains^^ 

Pip) = -Y,K-.\p)[^~w^^\p)K-.''{p)V- (19) 

As a result, knowing that F"(Q) = - K^a'^iQ), 
in the local approximation the irreducible polarizability 
P{q) for paramagnetic systems has the following familiar 
formii^ 

Pip) = p"ip)Hp) 

= P"(p)[l-|-t;,(q)g(p)P°(p)]-i (20) 

with the local-field factor Gip) = ip) / 2vcip) , where 

W'^-'^ip) = ^T^^W^-'^ip). This factor and the ex- 
change part of the local-field factor of Ref. 8 are formally 
the same. 

Next, we notice that, by representing the local in- 
teraction as W^'^ = Vc/e, the local-field factor can 
be expressed in terms of the RPA dielectric response 
function = 1 — VcP^ and the first order correction 
Ae(i) = i^c M%-'^ to 6° asM 

A similar expression for the imaginary part of Qiq) and 
with the longitudinal Lindhard dielectric function instead 
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FIG. 4: Feynman diagrams for the irreducible polarizability 
P in coordinate space. The RPA bubble (on the left) and the 
ladder diagrams (on the right) expressed in terms of the T 
matrix (shaded square) are represented here. 
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FIG. 5: Feynman diagrams for the direct (a) and exchange 
(b) terms of the electron self-energy. 

of £° was obtained in Ref. 0, where Ae^^^ contains the 
leading corrections to the RPA calculated within the 
model of the homogeneous electron gas. At w = 0, the 
factor H21|l is akin to the static local-field factor which 
has been calculated and parametrized in Ref. |2^ 

Thus, in the e—h case, we have the transparent connec- 
tion between the obtained local interaction and the ex- 
change part of the local-field factor arising from the first 
order in W exchange irreducible polarizability diagram. 
In this sense, the interaction W^~^ a gree s conceptually 
with the XC kernel considered in Ref. |3a. 



IV. SELF-ENERGY 

In this section we show how the electron self-energy 
and the T matrix 1)161) are related. As is known 
the electron self-energy obtained from the T matrix con- 
sists of a direct term and an exchange one (Feynman 
diagrams are shown in Fig. |SJ). The direct term 

E^(1,3)--*E J d2d4{G.,(4,2)T,V(l,2|3,4) 

+ G,,(2,4)r^V'(l,2|3,4)} (22) 

has e — e and e ~ h contributions, while the exchange 
term 

2^(2,3) =z/ dld4G,(4,l)T;-'=(l,2|3,4) (23) 
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FIG. 6: Four redundant diagrams originated from the T"^'" 
matrix (left column) and from the r'=~'' matrix (right col- 
umn). 



is defined by the spin-diagonal part of the T"^ matrix 
only. 

The Fourier transform of these terms leads to 



i^aa'i^^-.^^'P+k) 



, j.e-h( P+k P+k 

' (TO-' I r, I r, 



p^k)} (24) 



and 



Kip) 



i2ny 



dkGa{k) 

^, ye-e( k-p p-k 



P + k) (25) 



correspondingly. It is obvious from Eqs. H24() and H25() 
that with the T matrix of Eq. the exchange term 
and the spin-diagonal part of the e — e contribution in 
the direct term are, in fact, identical except for a sign. As 
a result, as well as in the Hubbard models, these terms 
are canceled. 

We notice here that, by substituting the T matrix as a 
solution of Eq. ifT^ into Eqs. if^ and if^ . one obtains^ 
for the direct term four lowest order diagrams (shown in 
Fig. O which disagree with the solution of the Hedin 
equationsi^ In order to avoid this problem, first of all, 
following Refs. |3 and we merely separate the first 
order exchange term (the GWA electron self-energy term 
j^Giv-j fj.Qj^ others. Next, we formally expand the T ma- 
trix into series, put into consideration a new value 
T"^, containing the second (or third in the e — h case) 
and higher order in W^^, items, and connect this value 
with the T matrix. This procedure yields 



(TO- 



T 

(7(7' 



(k) 

-\k) 



■pe — e 

: f' 



,J{k)Kl-^{k)Wp{k), 

-\k)[Ki-!^{k)w:-!^{k)f. 



<7<7' 



item in T^^^i 



On retaining the second order in W'^ 
we provide, thereby, the cancellation of the spin-diagonal 
e — e part of SlJ(p) and Yj%{p). Thus, additionally to the 
GW term, we obtain as a T-matrix contribution to the 
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FIG. 7: The imaginary part of the electron self-energy 
ImE[|q|, e(q)] of the electron gas as a function of momen- 
tum |q| at Ts ~ 2.07 (aluminium) and rs — 4.86 (potassium). 

ImE'"'^'^ andlmE" are shown by dashed and dotted lines, 
respectively. Solid line represents the sum of these terms. 
Insets: the electron lifetime t for the corresponding Vs val- 
ues, versus the excitation energy E — Ef- Dashed (solid) line 

shows r obtained from lmS°^ (Im[E°^ + E^"""]). e(q) is 
the free electron energy and qp (Ef) is the Fermi wave vector 
(energy). 



electron self-energy the following: 

S^b) = -(^ / dk{G-Ak)VZl{p + k) 

+ J2Ga'{k)T:-Mp-k)}. (26) 

Now we have only one term 

Kip) = 



dkG- 



r{k - p)w:-_iik) 

X K-_l{k)w^ 



which should be excluded from the T-matrix contribu- 
tion (|26|) . As a result, the electron self-energy can be 
expressed as So- = -I- — Yl^. The last item is 

an analog of the so-called double counting termiSiS4 In 
contrast to Ref. l2fil such item is present at the e — e 
contribution only. 

Employing the established connection between 
W^~^((l) and 5(g), one can, in principle, evaluate the 
ye-'i_niatrix contribution (|26|l (denoted as ) to 

the self-energy, additionally to the GW term, by using 
one of the local-field factors existing in the literature. 
But at present it can be seemingly done only for the 
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homogeneous electron gas for which these factors have 
been obtained and parametrized. 

Here, in order to roughly estimate the magnitude of 
YP"" , we exploit the static Q (q) of Ref. We have 
calculated the imaginary part of the electron self-energy 
for two values of the electron density corresponding to 
aluminium (r^ — 2.07) and potassium (vs — 4.86). Our 
results are shown in Fig. [3 It follows from the figure 
that in general S-^ is essentially less then S^iv espe- 
cially in the region where the decay due to creation of 
plasmons prevails. However, in the vicinity of the Fermi 
wave vector the T'^^'^-matrix contribution amounts on 
average to ~ 50% (70%) in relation to the GW term for 
Ts = 2.07 (4.86). This fact says that the contribution in 
question can be important in calculations of the decay 
of excited electrons whose initial energy is close to the 
Fermi energy. It is clear from the insets in Fig. |7| that 
the multiple scattering leads to shortening of the lifetime 
of such electrons. Note also that the values of the ratio 
/T,'-^'^ become greater when the electron density 
decreases. 



V. CONCLUSIONS 

In conclusion, we have presented a variational solution 
of the Bethe-Salpetcr equation which determines the T 
matrix describing multiple scattering both between two 



electrons or two holes and between an electron and a hole. 
The solution has been obtained within a local approxima- 
tion. The resulting expression for the T matrix is similar 
to that in Hubbard models but contains the local interac- 
tion depending on momentum and frequency. Thus the 
realized variational approach can be viewed as a method 
to evaluate the local interaction parameter U. In the case 
of multiple electron-hole scattering, a connection of this 
interaction with the local-field factors known from the 
literature has been established. We have also proposed 
a form of the T-matrix contribution to the electron self- 
energy which allows one to sum an infinite number of the 
electron-hole ladder diagrams for the electron self-energy 
without double counting. 
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